Impact of potentially toxic elements on pines in a former ore processing mine: Exploitation of hyperspectral response from needle and canopy scales.

Anthropic potentially toxic element (PTE) releases can lead to persistent pollution in soil. Monitoring PTEs by their detection and quantification on large scale is of great interest. The vegetation exposed to PTEs can exhibit a reduction of physiological activities, structural damage … Such vegetation trait changes impact the spectral signature in the reflective domain 0.4-2.5 µm. The objective of this study is to characterize the impact of PTEs on the spectral signature of two pine species (Aleppo and Stone pines) in the reflective domain and ensure their assessment. The study focuses on nine PTEs: As, Cr, Cu, Fe, Mn, Mo, Ni, Pb, Zn. The spectra are measured by an in-field spectrometer and an aerial hyperspectral instrument on a former ore processing site. They are completed by measurements related to vegetation traits at needle and tree scales (photosynthetic pigments, dry matter, morphometry …) to define the most sensitive vegetation parameter to each PTE in soil. A result of this study is that chlorophylls and carotenoids are the most correlated to PTE contents. Context-specific spectral indices are specified and used to assess metal contents in soil by regression. These new vegetation indices are compared at needle and canopy scales to literature indices. Most of the PTE contents are predicted at both scales with Pearson correlation scores between 0.6 and 0.9, depending on species and scale.

Mapping leaf metal content over industrial brownfields using airborne hyperspectral imaging and optimized vegetation indices.

Monitoring plant metal uptake is essential for assessing the ecological risks of contaminated sites. While traditional techniques used to achieve this are destructive, Visible Near-Infrared (VNIR) reflectance spectroscopy represents a good alternative to monitor pollution remotely. Based on previous work, this study proposes a methodology for mapping the content of several metals in leaves (Cr, Cu, Ni and Zn) under realistic field conditions and from airborne imaging. For this purpose, the reflectance of Rubus fruticosus L., a pioneer species of industrial brownfields, was linked to leaf metal contents using optimized normalized vegetation indices. High correlations were found between the vegetation indices exploiting pigment-related wavelengths and leaf metal contents (r ≤ - 0.76 for Cr, Cu and Ni, and r ≥ 0.87 for Zn). This allowed predicting the metal contents with good accuracy in the field and on the image, especially Cu and Zn (r ≥ 0.84 and RPD ≥ 2.06). The same indices were applied over the entire study site to map the metal contents at very high spatial resolution. This study demonstrates the potential of remote sensing for assessing metal uptake by plants, opening perspectives of application in risk assessment and phytoextraction monitoring in the context of trace metal pollution.

Unsupervised Monitoring Vegetation after the Closure of an Ore Processing Site with Multi-Temporal Optical Remote Sensing.

Ore processing is a source of soil heavy metal pollution. Vegetation traits (structural characteristics such as spatial cover and repartition; biochemical parameters-pigment and water contents, growth rate, phenological cycle…) and plant species identity are indirect and powerful indicators of residual contamination detection in soil. Multi-temporal multispectral satellite imagery, such as the Sentinel-2 time series, is an operational environment monitoring system widely used to access vegetation traits and ensure vegetation surveillance across large areas. For this purpose, methodology based on a multi-temporal fusion method at the feature level is applied to vegetation monitoring for several years from the closure and revegetation of an ore processing site. Features are defined by 26 spectral indices from the literature and seasonal and annual change detection maps are inferred. Three indices-CIred-edge (CIREDEDGE), IRECI (Inverted Red-Edge Chlorophyll Index) and PSRI (Plant Senescence Reflectance Index)-are particularly suitable for detecting changes spatially and temporally across the study area. The analysis is conducted separately for phyto-stabilized vegetation zones and natural vegetation zones. Global and specific changes are emphasized and explained by information provided by the site operator or meteorological conditions.

Monitoring oil contamination in vegetated areas with optical remote sensing: A comprehensive review.

The monitoring of soil contamination deriving from oil and gas industry remains difficult in vegetated areas. Over the last decade, optical remote sensing has proved helpful for this purpose. By tracking alterations in vegetation biochemistry through its optical properties, multi- and hyperspectral remote sensing allow detecting and quantifying crude oil and petroleum products leaked following accidental leakages or bad cessation practices. Recent advances in this field have led to the development of various methods that can be applied either in the field using portable spectroradiometers or at large scale on airborne and satellite images. Experiments carried out under controlled conditions have largely contributed to identifying the most important factors influencing the detection of oil (plant species, mixture composition, etc.). In a perspective of operational use, an important effort is still required to make optical remote sensing a reliable tool for oil and gas companies. The current methods used on imagery should extend their scope to a wide range of contexts and their application to upcoming satellite-embedded hyperspectral sensors should be considered in future studies.

Estimating persistent oil contamination in tropical region using vegetation indices and random forest regression.

The persistence of soil contamination after cessation of oil activities remains a major environmental issue in tropical regions. The assessment of the contamination is particularly difficult on vegetated sites, but promising advances in reflectance spectroscopy have recently emerged for this purpose. This study aimed to exploit vegetation reflectance for estimating low concentrations of Total Petroleum Hydrocarbons (TPH) in soils. A greenhouse experiment was carried out for 42 days on Cenchrus alopecuroides (L.) under realistic tropical conditions. The species was grown on oil-contaminated mud pit soils from industrial sites, with various concentrations of TPH. After 42 days, a significant decrease in plant growth and leaf chlorophyll and carotenoid contents was observed for plants exposed to 5-19 g kg-1 TPH in comparison to the controls (p < 0.05). Conversely, pigment contents were higher for plants exposed to 1 g kg-1 TPH (hormesis phenomenon). These modifications proportionally affected the reflectance of C. alopecuroides at leaf and plant scales, especially in the visible region around 550 and 700 nm. 33 vegetation indices were used for linking the biochemical and spectral responses of the species to oil using elastic net regressions. The established models indicated that chlorophylls a and b and ß-carotene were the main pigments involved in the modifications of reflectance (R2 > 0.7). The same indices also succeeded in estimating the concentrations of TPH using random forest regression, at leaf and plant scales (RMSE = 1.46 and 1.63 g kg-1 and RPD = 5.09 and 4.44, respectively). Four out of the 33 indices contributed the most to the models (>75%). This study opens up encouraging perspectives for monitoring the cessation of oil activities in tropical regions. Further researches will focus on the application of our approach at larger scale, on airborne and satellite imagery.

Application of PROSPECT for estimating total petroleum hydrocarbons in contaminated soils from leaf optical properties.

Recent advances in hyperspectral spectroscopy suggest making use of leaf optical properties for monitoring soil contamination in oil production regions by detecting pigment alterations induced by Total Petroleum Hydrocarbons (TPH). However, this provides no quantitative information about the level of contamination. To achieve this, we propose an approach based on the inversion of the PROSPECT model. 1620 leaves from five species were collected on a site contaminated by 16 to 77 g.kg-1 of TPH over a 14-month period. Their spectral signature was measured and used in PROSPECT model inversions to retrieve leaf biochemistry. The model performed well for simulating the spectral signatures (RMSE < 2%) and for estimating leaf pigment contents (RMSE ≤ 2.95 µg.cm-2 for chlorophylls). Four out of the five species exhibited alterations in pigment contents when exposed to TPH. A strong correlation was established between leaf chlorophyll content and soil TPH concentrations (R2 ≥ 0.74) for three of them, allowing accurate predictions of TPH (RMSE =3.20 g.kg-1 and RPD = 5.17). The accuracy of predictions varied by season and improved after the growing period. This study demonstrates the capacity of PROSPECT to estimate oil contamination and opens up promising perspectives for larger-scale applications.

Detection and discrimination of various oil-contaminated soils using vegetation reflectance.

The use of hyperspectral spectroscopy for oil detection recently sparked a growing interest for risk assessment over vegetated areas. In a perspective of image applications, we conducted a greenhouse experiment on a brownfield-established species, Rubus fruticosus L. (bramble), to evaluate the potential of vegetation reflectance to detect and discriminate among various oil-contaminated soils. The species was grown for 32 days on four different soils with mixtures of petroleum hydrocarbons and heavy metals. Additional plants were grown on either uncontaminated control or water-deficient soils for comparison. Repeated reflectance measurements indicated modified spectral signatures under both oil and water-deficit exposure, from leaf to multi-plant scales. The amplitude of the response varied with mixture composition, exposure time, acquisition scale and spectrum region. Reflectance changes were linked to alterations in chlorophyll, carotenoid and water contents using vegetation indices. These indices were used to catch spectral similarities among acquisition scales and to discriminate among treatments using Kendall's coefficient of concordance (W) and regularized logistic regression. Of the 33 vegetation indices tested, 14 were concordant from leaf to multi-plant scales (W > 0.75, p < 0.05) and strongly related to leaf biochemistry (R2 > 0.7). The 14 indices allowed discriminating between each mixture and the control treatment with no or minor confusions (≤5%) at all acquisition scales, depending on exposure time. Some of the mixtures remained difficult to discriminate among them and from the water-deficit treatment. The approach was tested at the canopy scale under natural conditions and performed well for identifying bramble exposed to either one of the experimentally-tested mixtures (90% accuracy) or to uncontaminated soil (83% accuracy). This study provided better understanding of vegetation spectral response to oil mixtures and opens up promising perspectives for future applications.

Assessing Soil Contamination Due to Oil and Gas Production Using Vegetation Hyperspectral Reflectance.

The remote assessment of soil contamination remains difficult in vegetated areas. Recent advances in hyperspectral spectroscopy suggest making use of plant reflectance to monitor oil and gas leakage from industrial facilities. However, knowledge about plant response to oil contamination is still limited, so only very few imaging applications are possible at this stage. We therefore conducted a greenhouse experiment on three species long-term exposed to either oil-contaminated or water-deficient soils. Reflectance measurements were regularly performed at leaf and plant scale over 61 days of exposure. Results showed an increase of reflectance in the visible (VIS), the red-edge and the short-wave infrared (SWIR) under both oil and water-deficit stress exposure. A contrasted response in the near-infrared (NIR) was also observed among species. Spectra underwent transformations to discriminate species' responses to the different treatments using linear discriminant analysis (LDA) with a stepwise procedure. Original and transformed spectra enabled to discriminate the plants' responses to the different treatments without confusion after 61 days. The discriminating wavelengths were consistent with the spectral differences observed. These results suggest differential changes in plant pigments, structure and water content as a response to various stressors, and open up promising perspectives for airborne and satellite applications.

Estimation of soil moisture content from the spectral reflectance of bare soils in the 0.4-2.5 µm domain.

This work aims to compare the performance of new methods to estimate the Soil Moisture Content (SMC) of bare soils from their spectral signatures in the reflective domain (0.4-2.5 µm) in comparison with widely used spectral indices like Normalized Soil Moisture Index (NSMI) and Water Index SOIL (WISOIL). Indeed, these reference spectral indices use wavelengths located in the water vapour absorption bands and their performance are thus very sensitive to the quality of the atmospheric compensation. To reduce these limitations, two new spectral indices are proposed which wavelengths are defined using the determination matrix tool by taking into account the atmospheric transmission: Normalized Index of Nswir domain for Smc estimatiOn from Linear correlation (NINSOL) and Normalized Index of Nswir domain for Smc estimatiOn from Non linear correlation (NINSON). These spectral indices are completed by two new methods based on the global shape of the soil spectral signatures. These methods are the Inverse Soil semi-Empirical Reflectance model (ISER), using the inversion of an existing empirical soil model simulating the soil spectral reflectance according to soil moisture content for a given soil class, and the convex envelope model, linking the area between the envelope and the spectral signature to the SMC. All these methods are compared using a reference database built with 32 soil samples and composed of 190 spectral signatures with five or six soil moisture contents. Half of the database is used for the calibration stage and the remaining to evaluate the performance of the SMC estimation methods. The results show that the four new methods lead to similar or better performance than the one obtained by the reference indices. The RMSE is ranging from 3.8% to 6.2% and the coefficient of determination R2 varies between 0.74 and 0.91 with the best performance obtained with the ISER model. In a second step, simulated spectral radiances at the sensor level are used to analyse the sensitivity of these methods to the sensor spectral resolution and the water vapour content knowledge. The spectral signatures of the database are then used to simulate the signal at the top of atmosphere with a radiative transfer model and to compute the integrated incident signal representing the spectral radiance measurements of the HYMAP airborne hyperspectral instrument. The sensor radiances are then corrected from the atmosphere by an atmospheric compensation tool to retrieve the surface reflectances. The SMC estimation methods are then applied on the retrieve spectral reflectances. The adaptation of the spectral index wavelengths to the HyMap sensor spectral bands and the application of the convex envelope and ISER models to boarder spectral bands lead to an error on the SMC estimation. The best performance is then obtained with the ISER model (RMSE of 2.9% and R2 of 0.96) while the four other methods lead to quite similar RMSE (from 6.4% to 7.8%) and R² (between 0.79 and 0.83) values. In the atmosphere compensation processing, an error on the water vapour content is introduced. The most robust methods to water vapour content variations are WISOIL, NINSON, NINSOL and ISER model. The convex envelope model and NSMI index require an accurate estimation of the water vapour content in the atmosphere.

Aggregation process of optical properties and temperature over heterogeneous surfaces in infrared domain.

We propose a modeling of the aggregation processes of optical properties and temperature over the heterogeneous landscape in the infrared domain (3-14 microm). The main objectives of the modeling are to understand how these parameters aggregate and to study their links at different spatial scales. As the landscape is described at each scale by its radiative parameters, general equations linking the radiative parameters at a given high spatial scale to those at a rough scale are proposed. Then these equations are applied to several synthetic landscapes. An analysis based on a design of experiments is conducted to point out the influence of each of the input factors. The results show the importance of the intrinsic parameters (reflectance, emissivity, and surface temperature) of each surface element and also the directional and spectral behaviors of the aggregated parameters.

Thermal infrared radiance simulation with aggregation modeling (TITAN): an infrared radiative transfer model for heterogeneous three-dimensional surface--application over urban areas.

The thermal infrared radiance simulation with aggregation modeling (TITAN) model, presented here, is an innovative transfer radiative code in the infrared domain (3-14 microm). It takes into account the three-dimensional (3D) structure of the landscape and simulates all the radiative components introduced by this 3D structure, which are due to the reflection and emission of walls and sloping roofs. Examples are given to illustrate the new opportunities offered by TITAN over urban areas. First, a phenomenological study is conducted at four wavelengths analyzing the relative effect of all the radiative contributors to the total signal. The same analysis is performed at bottom of atmosphere, which reveals an error occurring when a flat assumption is made (between 1% and 5%). In a second example, the directional effects at sensor level are simulated and show that the radiative temperature can vary by up to 10 K.